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Abstract 

Some practical improvements are proposed for the "optical-shaker" laser-cooling technique 
[I. S. Averbukh and Y. Prior, Phys. Rev. Lett. 94, 153002 (2005)]. The improved technique 
results in an increased cooling rate and decreases the minimum cooling temperature achievable 
with the optical shaker. The modified shaker requires only one measurement of the force on the 
atoms before each cooling step, resulting in a simplification of the feedback electronics. The force 
is inferred from the power variations of the transmitted laser beams and is used to determine the 
best moment at which the cooling steps are applied. The temperature of the atomic system is 
automatically monitored, which allows maintaining an optimum cooling rate as the temperature 
decreases. The improved shaker is simple to build, provides a faster rate of cooling, and can work in 
the microkelvin regime. Numerical modeling shows a reduction by a factor of three in the required 
number of phase jumps and a lower temperature limit reduced by an order of magnitude compared 
to the initially proposed shaker. The technique is also extendable to cooling in three dimensions. 
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I. INTRODUCTION 



Laser cooling of atoms [l[ has made possible important experiments such as high-precision 
measurements of fundamental constants and the observation of Bose-Einstein condensation. 
In many experiments, laser cooling relies on the nearly resonant exchange of momentum 
between the optical field and the atoms to obtain significant cooling. However, conventional 
laser cooling is restricted to atomic species which offer an optical cycling transition as a 
nonconservative force. Nonresonant interactions can address a larger variety of atoms and 
molecules by providing the means of trapping these particles and isolate them from strong 
resonant interactions with the lasers. One disadvantage is that the conservative nature of 
the force does not naturally lend itself to cooling through an energy loss mechanism. Only a 
few cooling techniques usingnonresonant laser interactions have been proposed 
Optical lattice clocks jf], Q, 3, are good examples where weakly interacting laser beams 
isolate the trapped atoms from external perturbations. Quantum computers also require 



cold atoms isolated from external influences to minimize decoherence 
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One example of the nonresonant cooling technique is the optical shaker proposed by 



Averbukh and Prior 
and Sisyphus cooling 



(AP shaker). It combines the ideas of stochastic laser cooling 
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13( | without the need for a resonant laser interaction. Both mecha- 



nisms contribute to remove energy from trapped particles. Stochastic cooling relies on the 
statistical nature of the spatial distribution of the particles in an optical potential to reduce 
their potential energy. Such an optical potential can be produced with two counterprop- 
agating laser beams. For particles in a sinusoidal potential, the total potential energy of 
the particles E is nonzero due to the random spatial distribution of the particles. This 
energy E is a sinusoidal function of the position of the optical potential. By imposing a 
rapid phase jump on one laser beam at the appropriate time, the optical potential can be 
spatially displaced to decrease the potential energy of the particles without changing their 
kinetic energy. The second mechanism, Sisyphus cooling, decreases the total kinetic energy 
of particles when they climb the potential hills. The potential energy E increases with time 
until another phase jump decreases it. Repeated applications of phase jumps lowers the total 
energy of the system and reduces its temperature. In stochastic cooling, the particles are not 
thermalized by collisions. Throughout this paper, collisions are ignored and an equivalent 
temperature proportional to the total kinetic energy is used. 
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The optical shaker produces phase jumps at the appropriate time and with the appro- 
priate amplitude. For the atoms to lose kinetic energy, a sufficient amount of time is given 
to let the system evolve so that the particles will have climbed the potential hills before 
the next phase jump. A fixed time interval is chosen based on the wavelength of the opti- 
cal potential and the average velocity of the particles. The required size of the phase step 
is more difficult to evaluate. As is explained below, the position of the optical potential 
must be inferred from the difference in power of the transmitted laser beams. Before each 
jump, two consecutive measurements are required from which the size of the phase step is 
calculated numerically. Fast electronics and a fast processor are required since under typical 
experimental conditions, only a few microseconds are available for the measurements and 
the calculation of a trigonometric function. One drawback of the method is that it is limited 
to a minimum temperature in the several tens of microkelvins and the time constant to reach 
this temperature is of the order of several tens of seconds. At this temperature, the power 
signals become very small and the cooling rate is limited by the shot noise on the detectors. 

The modified optical shaker proposed in the present paper leads to simplifications of the 
experimental system, an increase of the efficiency of the cooling process, and a lower cooling 
temperature limit. To achieve this, the modified shaker monitors the time dependence of 
the transmitted powers. The modified shaker is more efficient because it directly monitors 
the force on the particles instead of the potential energy. A phase jump is triggered at 
the most appropriate time to achieve the fastest cooling. Only one measurement is required 
before each jump, therefore reducing the required bandwidth of the electronic circuits. Since 
the size of the phase step is a linear function of the time between the jumps, the need for 
a complex processor is eliminated. With these improvements, atoms can be cooled down 
to temperatures near one microkelvin. The theoretical basis describing the operation of 
the modified shaker is essentially the same as for the originally proposed AP shaker J4]. 
This theory is presented in Sec. Hi] with a calculation of the cooling rate under optimized 
conditions. The modified shaker is described in Sec. II III with a theoretical derivation of it's 
efficiency. In Sec. IIV} both a numerical and an analytical derivation of the temperature as 
a function of the cooling time are presented. A comparison is made between the numerical 
model and the analytical model. The limitations resulting from shot noise on the detectors 
and the minimum achievable temperature are discussed. 
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II. FEEDBACK MINIMIZING THE POTENTIAL ENERGY 

A qualitative description of the AP shaker [4] can be given as follows. Consider first the 
simple model of a single particle moving along two counterpropagating laser beams. A dipole 
force on the particle arises from the momentum kick when a photon is taken from one beam 
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151 ]. Since the force depends 



and transferred to the other beam via stimulated emission 
on a transfer of light from one beam to the other, it is monitored directly by measuring the 
power difference between the two beams. The use of a periodic sinusoidal potential in the 
description of the mechanism is justified since the force is conservative. To cool the particle, 
its kinetic energy must be reduced. This can be achieved in two steps. First, the potential 
energy of the particle is reduced with a rapid translation of the optical potential to a position 
where the particle will be at a minimum. The second step consists in waiting a sufficient time 
to let the particle move away from a minimum and lose kinetic energy. These steps can be 
repeated numerous times to reduce considerably the kinetic energy of the particle. Consider 
now N particles moving along the optical potential. Because the particles have random 
positions, the total potential energy of these particles will not completely vanish. The total 
potential energy will be sinusoidally dependent on the position of the optical potential. As 
in the case of one particle, one can rapidly translate the optical potential to a position where 
the total potential energy is minimized. The system will then evolve and the particles, as 
an ensemble, will lose some kinetic energy and therefore reduce their temperature. The 
velocity is assumed to be a thermal distribution, but the spatial distribution of the particles 
is nearly uniform across one wavelength of the periodic potential. The spatial distribution 
is essentially independent of time, except immediately after a phase jump when there is a 
larger density near the potential minima. 

Quantitatively, the AP shaker operates on an ensemble of N particles interacting with 
two linearly polarized and counterpropagating laser beams having an electric field amplitude 
Eocos(kiz — uit + <p) and Eq cos (kiz + Uit) , where ki = 2n/\. The beams have the same 
polarization and the same frequency uj\ = 2-kc/X with a phase difference <fi at z = 0. The 
total potential energy of the N particles at positions Zi(t) is written 

U(4>,t) = -U u(4>,t)/2, (1) 

where the normalized potential energy is u((p, t) = X)i=i cos[zi(t)+(f)), the normalized position 
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is Zi(t) = 2irZi(t) / A, and A = A/2. The total force on these atoms is then 

F^,t) = -nU f(<P,t)/A, (2) 

where the normalized force /(</>, t) = Y^iLi sia[zi(t) + 0], Uq = aE^, and a is the nonresonant 
polarizability. In what follows, a red-detuned laser beam is considered so that a > 0. 
Stochastic cooling uses the random fluctuations in the potential energy of the system it(0, t) 
to decrease the temperature of the system. By imposing a rapid phase jump from to 
+ 50, the total potential energy is decreased without changing the kinetic energy of the 
particles. For a given it(0, t) and /(0, t), the required phase step to bring the potential 
energy U((f) + 50, t) to a minimum is [4] 

5(f) = arctan (— /(0, t), +w(0, t)) , (3) 

where the notation arctan (psin(#), pcos(6)) = 9 (for any p > and — 7r < < 7r) is used to 
resolve the angle ambiguity encountered with the regular arctan (sin(#)/ cos(#)) function. A 
phase jump causes the system to lose an energy equal to 

U 



NAE 

2 



+ v / M 2 (0,t) + / 2 (0,t) - u(<f>, t) 



(4) 



where the positive root is used 



4| . The efficiency of the cooling technique is quantified by the 
average energy loss per particle at each jump, (AE) . Since the system is allowed to evolve for 
a sufficiently long time before a phase jump, the steady-state statistical distributions of the 
two variables u(<p,t) and /(0, t) are used for the calculation of (AE) . Random fluctuations 
cause these two variables to have a normal distribution with a standard deviation 



a = y/N/2. (5) 



Their average values are more difficult to evaluate because the spatial distribution of the 
particles is not a uniform random spatial distribution. Fast particles tend to accumulate 
near the potential maxima where they are slower, while slow particles are trapped at the 
bottom of the potential wells. The result is a sinusoidally modulated distribution in phase 
with the optical potential. In general, the expectation value (U((j),t)) is not zero. However, 
since the force is in quadrature with the potential, it is also in quadrature with the spatial 
modulation of the particle distribution and averages to a null value 

(F(0,t)) = O. (6) 



To evaluate the average potential {U(<f>, t)) , the spatial distribution of the particles needs 
to be known. Slower particles near the top of the potential hills spend more time in a region 
of high potential than faster particles near the bottom of the potential wells. The time- 
dependent potential energy is therefore positive for longer time periods than it is negative. 
As a result, the average potential energy {4] is 

(U(<P,t)) s = rNU /2 (7) 

for Uq << ksT, where r = Uq/ '(4/cbT) is a parameter proportional to the cooling laser 
intensity. This relation is valid for a specific velocity class and also for a Maxwell-Boltzmann 
speed distribution (hence the subscript S). However, when the v z component of a Maxwell- 
Boltzmann velocity distribution is considered, one finds that a large fraction of the particles 
have a small velocity component v z and remain trapped in the potential wells. These slow 
particles have a negative time-averaged potential energy if their total energy (kinetic and 
potential) is smaller than = (0.327. . -)Uq. Particles with a total energy above Ef. have 
a positive average potential energy, but that energy is not high enough to compensate 
the negative contribution from the larger number of slower particles. The time-averaged 
potential energy can be estimated for a thermal distribution of the velocity components 
in equilibrium with the standing-wave potential. This requires an integral done over the 
potential energy e(z) = —UqCos(z)/2 of a particle at location z multiplied by the energy 
distribution function of the number of particles at that location: 

The average potential energy for a Maxwell-Boltzmann velocity distribution is 

(U(4>,t)) v = -rNUo/2 (8) 

for a small potential Uo such that r << l/\ / 2N. Note the negative value of the averaged 
potential compared to Eq. (j7]). This negative potential energy is a result of the maximum 
value at v z = of the Maxwell-Boltzmann velocity distribution. For r » 1 /\ / 2N, the 
average potential energy reaches the limit (U((j),t)) = —NUq/2. This corresponds to having 
all the particles at the bottom of the potential wells. Averaging the energy loss given by 
Eq. (JI]) over the statistical distributions and the average values of /(</>, i) and u(<f),t) [from 
Eqs. (JSJ), (jUJ) and (jSJ)], the average energy loss per particle at each jump [4J is 

(AE) = k B T/(2N) (9) 
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FIG. 1: Schematic setup of the modified shaker. The shaker beams are indicated by the long 
black arrows. If separate trapping beams are used, they must be configured in a way that will not 
interfere with the shaker beams (gray arrows) . Two detectors measure the transmitted powers from 
which the power difference AP is calculated. A zero-crossing detection circuit triggers the phase 
jump with an amplitude calculated from the slope of the derivative of AP and the time interval 
since the last phase jump. The output of a rate meter controls the laser intensity for optimum 
cooling speed. 



and is independent of U as long as the condition r > l/\/2N is satisfied. This result 

assumes a Maxwell-Boltzmann velocity distribution of the particles. Collisions are ignored 
in this model. 

Experimentally, cold particles are trapped by two nearly counterpropagating laser beams 
as shown in Fig. [TJ The standing wave provides the potential wells used for the stochastic 
cooling. An electro-optic modulator (EOM) can rapidly change the phase of one of the 
beams. Two detectors measure the transmitted power of the two beams from which the force 
on the atoms f(<p, t) is calculated by monitoring the difference. Conservation of momentum 
imposes that AP(</>, t) = Pi — P 2 = cF(0, t) = —TrcU f(4>, t)/A, where c is the speed of light 
and Pi and P2 are the powers measured by the detectors. An estimate of the stochastic 
power exchange using Eq. ([5]) gives 



The AP shaker also needs the measurement of u(<f),t) to calculate the size of the necessary 
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phase step with Eq. ([3]). The potential energy is measured by briefly imposing a 7r/2 phase 
step and measuring AP(cf) + n/2,t) = cF(<fi + n/2,t) = — 7icUou((fi, t)/A. The phase step 
calculated from Eq. ([3]) is then added to the phase <fi to increase the number of particles 
near the bottom of the potential wells. According to the mechanism of Sisyphus cooling, 
the particles will climb the potential hills, gain potential energy, and lose kinetic energy. 
Successive applications of the appropriate phase steps 5(f) will reduce the total energy of the 
system, therefore resulting in cooler particles. 

In order to give enough time for the system to evolve sufficiently before a phase jump, the 
time between jumps tj is chosen to be near the time of an oscillation period of low-energy 
particles trapped by the wave tn = (2r)~ 1 ' 2 t\, where t& = A/v rms is a characteristic time 



of the evolution of the system and v rms = ^hsT/m is the root-mean-square velocity of 
the particles at temperature T in one dimension. The laser intensities are chosen so that 
r ~ 1/ \ / 2N. This choice of parameters will result in an average energy loss given by Eq. (J9]) 
and also produce reasonably sized phase steps <50 of the order of one radian. A simpler 
feedback scheme j^j] is possible where only a single measurement of the sign of AP(<f), t) is 
needed. A fixed phase step ±<50o is applied to displace the standing wave opposite to the 
direction of the instantaneous force. In this case, however, a phase jump will result in a 
smaller average energy loss per particle of {AE) = ksT '/(ttN). Adaptive lowering of Uq is 
indicated to keep the optimum cooling conditions. 

The results of the example given in 4j with r = 0.05 and tj = tn are reproduced in 
Fig. 12(a), showing the temperature as a function of the number of phase jumps. It is also 
interesting to know how quickly the cooling occurs by plotting the temperature as a function 
of time in Fig. [3](a). An initial temperature of 100 fiK is used in this example. One can 
ask how much time is needed between jumps in order to achieve the fastest cooling rate 
- that is, to maximize the energy loss per unit time and yet provide enough time for the 
system to evolve and mix sufficiently. One could argue that some time after a jump, the 
system reaches a state of highest potential energy and lowest kinetic energy when most of 
the particles go from the bottom to the top of a potential hill. The time required for this 
is about tj = t\/2, which is much shorter than £q. Numerical simulations show that the 
efficiency of the phase jumps remains constant for times as short as tj pa 0.22t\, but that 
the cooling rate is improved by a large factor. As suggested in {4], a decreasing laser power 
allows a more efficient cooling as the temperature decreases. The numerical simulation was 
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FIG. 2: Temperature as a function of the number of phase jumps. In (a), the phase step of the AP 
shaker is determined from two measurements. Curve (b) is the same technique with adaptive laser 
power and a shorter time between jumps. Curve (c) shows the cooling using the method proposed 
in this paper. The curves are the result of an average of fifteen runs. 

repeated with a decreasing laser power kept at a level proportional to the temperature and 
a faster jump rate tj. Optimization of the power gives the fastest cooling rate with an initial 
r « 0.07 for N = 1000 particles. The results are plotted in Figs. 12(b) and [3(b), showing 
a net improvement of the cooling efficiency as a function of time. With a varying laser 
intensity, the cooling process follows Eq. (Q over a wide range of temperatures. 

III. FEEDBACK MINIMIZING THE FORCE 

A simplification of the experimental setup resulting in a faster cooling rate is presented 
in this section. To understand how the improved technique works, consider again the simple 
system with a single particle rolling up and down the potential hills. When the particle 
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reaches the top of the hill, the largest energy loss is obtained with a phase step 5(f) = tt, 
which brings the particle to the bottom of a well. As discussed above, it is difficult to know 
when the particle reaches that position since one can only directly measure the force on 
that particle F(<ft,t), not the potential U(<p,t). A possible alternative is to wait until the 
particle reaches half the height of the potential well, then impose a phase step 50 = ±7r/2. 
This has the advantage that the force reaches a measurable extremum when it is time for 
a phase jump. The sign of the phase step is readily determined from the sign of the force. 
The process is different when there are several particles in the system since the potential 
energy does not oscillate in quadrature with the force, but instead is an independent random 
process. However, the technique is still applicable with a small modification. 

For reasonably small phase steps \5<p\ < 1 near a potential minimum, the relation 
\f((f>,t)\ « \u(<p,t)\ is valid. From Eqs. ([I]) and ([8]) we get (u((p,t)) v = rN. This im- 
plies that u(<p,t) is often positive. In this case, Eq. (J4]) simplifies to 



When a phase jump is applied, the largest energy loss is obtained with the largest possible 
value of f 2 (<ft,t) as in the one-particle example given above. By monitoring the force and 
waiting until it reaches an extremum, the phase jump actually produces the largest energy 
loss. Replacing u((p,t) by its average value rN and replacing f 2 ((p,t) by the variance N/2, 
Eq. (TIT]) gives (AE) ~ k^T / '(2N), the same energy loss obtained with the more complex 
method using measurements of the potential energy and the force on the particles given 
by Eq. ([9]). Because only one measurement is needed instead of two, the bandwidth of the 
detection system can be halved, resulting in an improved signal-to-noise ratio. Additionally, 
the time between jumps is automatically controlled by the measurement of the force and 
automatically adapts itself to the dynamics of the system. 

One must determine the sign and magnitude of the phase step that will provide the best 
cooling. After some time At the force reaches an extremum and assuming that u(4>,t) is 
positive, the sign of the phase step is inferred from the sign of f(4>,t) and Eq. ([3]). Since 
no direct measurement of the potential u{cj),t) is available, the magnitude of the required 
phase step cannot be calculated. One could use a constant value for the phase step 5(p, 
but this would reduce the efficiency of the cooling process and is not adaptable to different 
temperatures. One could also assume a constant value of the potential, but an accurate 
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FIG. 3: Temperature as a function of time. In (a), the phase step for the AP shaker is determined 
from two measurements. The time between phase jumps is equal to the period of oscillation of 
low-energy particles trapped by the wave. Curve (b) is the same technique with adaptive laser 
intensity optimized for the fastest cooling rate and a shorter time between jumps tj ps 0.22t\. In 
(c), cooling is obtained using the modified shaker. Because there is no need to let the system freely 
evolve and remix between two consecutive phase jumps, the cooling time is significantly shorter. 
The temperature follows the expected theoretical curve shown as a thin line after 9 ms. The curves 
are the result of an average of fifteen runs. 

measurement of f(<f>,t) and a calculation of the arctangent is required. An easier solution 
is to use the measured time At elapsed since the last phase jump. Just after a phase jump 
which brings the force to f((j),t) rs 0, the important point to notice is that the force evolves 
according to a random walk, which causes the value of f 2 ((p, t) to increase linearly with time. 
The required phase step also increases linearly with time as 



5(j) = ±u(r,v rms )AT. 



(12) 
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The equations of motion show that the time evolution of the system is inversely proportional 
to the velocity v rms . This suggests the use of a variable phase step factor u(r, v rms ) = 
<t>{T)v rms / h, where the proportionality factor <f>(r) is not explicitly dependent on v rms . 



IV. NUMERICAL RESULTS AND DISCUSSION 



Using numerical simulations, the fastest cooling is obtained with the optimized param- 
eters r 0.098 and <p(r) = 0.85 for N = 1000 particles. The good performance of the 
technique is readily visible in Figs. [2(c) and[3](c). A 40% reduction in cooling time is ob- 
tained compared to the optimum adaptive cooling rate obtained with the AP shaker (4]. 
With optimum cooling conditions, the average size of a phase step {\8<f>\) ~ 0.21 satisfies the 
above requirement of a small angle. Since the jumps occur when f 2 (<f),t) is at a maximum, 
a larger AE per jump is obtained. The average energy loss per particle per jump becomes 
as large as 

(AE) « lAk B T/(2N). (13) 

The 40% improvement over the estimated energy loss results from the fact that the variance 
of f(4>,t) is actually larger than A^/2 at the time of the jump. The larger variance shows 
that the system has not yet evolved to a random distribution just before the jump. One 
sees the net advantage in choosing the appropriate time before imposing the phase jump 
as opposed to evenly distributing the jumps in time. With a larger number of atoms, the 
optimum energy loss per jump given by Eq. (TT~3]) is achieved by scaling down the power of 
the shaker beams as \/yN. This maintains the average size of the phase steps near the 
optimum value of 0.21 rad. 

Experimentally, the technique can be implemented in the following way. In order to 
detect the extremum of the force on the atoms, a differentiating circuit is used to process 
the signal AP. This is indicated by d/dt in Fig. [TJ The circuit is easily built with a high- 
frequency pass RC filter. The extremum occurs when dAP/dt = and is measured with 
a zero-crossing detection circuit. Because of the increased high-frequency noise that could 
result at the output of the filter, a high-frequency cutoff filter should be added. Numerical 
simulations showed that a filter with a response time t& is a good choice. Three signals are 
then derived from the zero-crossing detection circuit: (a) a signal proportional to the time 
interval At since the last zero crossing (i.e., a ramp generator), (b) the sign of the derivative 



12 



just after the zero crossing (i.e., a comparator), and (c) a trigger signal for a rate meter. 
When multiplied together, the first two signals give the value of 5(f) to be added to the phase 
of the EOM as per Eq. (1121) . The rate meter measures the average number of phase jumps 
per unit time interval and gives a signal proportional to the velocity u\ = Cnv rms . Numerical 
simulations show that Cr = 4.2/A. 

The equations of motion also show that to keep the system optimized, the parameter 
r = Uo/ \AksT) should be kept constant at the value determined numerically as r « 3.1/ y/N. 
This imposes that for a given number of atoms, the potential Uo remains proportional to 
the temperature as a function of time - that is, proportional to v\. This value can easily be 
obtained with a squaring circuit connected to the output of the rate meter. This technique 
can easily be implemented as an all- analog circuit, therefore making it easy to achieve a fast 
response time of 100 ns needed for initial temperatures of the order of 10 mK. It also has 
a net advantage over the method requiring two measurements and rapid digital processing 



for each phase jump |4|, [17(|. With this adaptive response implemented, the rate of the 
phase jumps is proportional to v rms or \/T and the energy loss follows Eq. ( fl3l) . This is 
confirmed by the results plotted in Fig. G](c), where the thin line extended after 9 ms shows 
the temperature of the system calculated with the differential equation 

dT(t)/dt « -l.AC R v rms T{t)/N oc -T(tf 2 /N. 

The improved performance of the modified shaker is demonstrated here with a few exam- 
ples. A cooling time of the order of several seconds is reported for the original AP shaker 4| 
with the parameters N = 10 6 cesium atoms initially at 15 mK, a laser power P = 2 fiW 
focused to a beam waist wq = 10 /im and A = 1 fim with a detuning A = 100r, where T is 
the decay rate of the upper level of a two-level atom. The cooling time is reduced to 300 ms 
with the modified shaker method. (In this case, the optimum power is 0.12 fiW at 15 mK, 
and the cooling time is defined as the time required to reach e~ 2 ~ 13.5% of the initial 
temperature.) Another advantage of the improved technique is that it does not depend on 
an accurate measurement of the dc power of the transmitted laser beams (the AP shaker re- 
quires a power balance good to a few percent). Instead, the modified shaker measures power 
fluctuations. This is useful since at colder temperatures the stochastic power exchange AP 
becomes small and difficult to measure [3] as it scales proportionally to the intensity of 
the laser beams, which is adaptively decreased to maintain optimum cooling. The use of a 
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time differentiating circuit in the modified shaker improves performance since the output is 
insensitive to a possible dc signal which could appear, for example, as a result of interference 
from reflections on surfaces. 

The lowest temperature achievable with this technique is limited by the detector's shot 
noise and the intensity noise on the laser beams. The effects of intensity noise have been 
discussed previously 17|. The shot-noise limit for the modified shaker is reduced because the 
detectors operate with half the bandwidth of the AP shaker. Thus, colder temperatures can 
be achieved. The signal to be detected can be estimated 4J as AP/P ~ 4y / 2iV oc / (Aeoifo) , 
where a ~ 47reo(c/co>o) 3 (r/A). For the best cooling efficiency, the conditions derived above 
require a shaker power P oc WqT(A/T)/\^N. This results in a detected signal AP oc T. 



On the other hand, the shot-noise power is Pn oc y/PrjB, where B oc vT is the required 
bandwidth and rj is the number of measurements needed for the shaker algorithm. At the 
minimum cooling temperature, these two powers are equal and 

w^ 2 /A x 2 
T min « 0.018 fjK 1)1 ' 



n vr 

where Wo is the beam waist in fim and the constant is found from the optimum parameters 
calculated above. In the example above for the AP shaker, 77 = 2 and the power difference 
becomes equal to the quantum noise at T m i n = 41 /zK. The modified shaker requires only 
half of the bandwidth used by the AP shaker. With rj = 1 the minimum temperature is 
reduced to T min = 10 /iK. The time required to cool from 74 fiK to 10 /iK is 4 s. These 
results are consistent with the conclusions that the minimum temperature obtainable with 
an optical shaker is not limited by the photon recoil limit, but by the square of the resolution 
of the measurement [ijij]. If the coldest temperature is the goal to be achieved, a smaller 
detuning and a small beam waist should be considered. 

As the atoms are cooled, the optimum cooling power will decrease and the trapping 
potential will become shallower. If the shaker beams are used as trapping beams, the trap 
depth might not be sufficient to hold the atoms against gravity. However, the shaker beams 
can be distinct from the trapping beams. The operation of the optical shaker is possible 
with separate trapping beams, as long as the standing waves produced by the beam pairs 
do not spatially coincide and that no power from the trapping beams reaches the power 
detectors. The first condition is accomplished by using an angle between the two types 
of beams. The second condition is achieved with a frequency offset to avoid producing a 
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standing wave between the two types of beams. The power in the shaker beams can be 
decreased independently of the power in the trapping beams, therefore enabling efficient 
cooling. This configuration is depicted in Fig. [TJ It also has the significant advantage that 
a lower cooling temperature is achievable. If a small angle is used between the shaking 
and tapping beams, the atoms see a potential with an effective wavelength A longer by a 
factor esc (8) in a direction perpendicular to the trapping beams. The atoms have a larger 
amplitude of oscillation in that direction, which has the effect of decreasing the value of Cr. 
The cooling speed is slower, but since the phase jump rate is decreased, a smaller bandwidth 
can be used for the detectors, therefore decreasing the minimum achievable temperature. 
For example, an angle 6 ~ 0.32 rs 18.4° decreases the cooling rate by a factor of vTo, but 
reduces T min by a factor of ten. In some applications the trapping beams have to be tuned 
very far away from any resonant level. However, the cooling beams can be detuned by only 
a few GHz for the optimal cooling efficiency and the coldest temperature. 

The technique is easily generalized to an n-dimensional shaker built using a minimum of 

[ I 

n + 1 shaker beams [20J. In that case, phase control of n beams is required based on n + 1 
power measurements. In general, the phase control on beam number i of an n-dimensional 
shaker uses the derivative of APj = (X^yi-P?) — Pi- If trapping beams are present, they 
must allow some motion of the particles in at least one dimension. 



V. CONCLUSIONS 



A practical and efficient technique to implement optical shaking for laser cooling is de- 
scribed. The technique combines the simplicity of the electronic feedback system and the 
adaptability required for efficient cooling over a large temperature range. The faster re- 
sponse time of the feedback loop improves the limit imposed by its finite response 1?J • The 
modified shaker does not let the system go back to a random distribution between each 
shake, therefore improving the cooling speed compared to the AP shaker 4j. The reason 
behind the improvement is that this technique is based on minimizing the force instead of 
the potential energy, which takes advantages of the correlation in the distribution of the 
atoms immediately after a cooling step. By responding to immediate changes in the power 
fluctuations, it provides a high cooling rate that exceeds the rate obtained when predeter- 
mined time intervals are used between jumps. This makes the operation of the improved 
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shaker possible down to tens of microkelvins. When used with additional trapping beams, 
lower temperatures can be achieved. A fast cooling rate can become especially important 
in an optical lattice clock where the atoms need to be kept in the trapping potential at the 
coldest possible temperature for long periods of time. The simplicity of this scheme makes 
it a promising configuration for a demonstration of the modified shaker cooling. 
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